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Abstract 

A representation of divisor function r(n) = o"o(n) by means of logarithmic residue of a 
function of complex variable is suggested. This representation may be useful theoretical 
instrument for further investigations of properties of natural numbers. 
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1 Introduction 



Properties of natural numbers are the fundamental question of mathematics. Investi- 
gations of these properties in the Ancient Times (e.g. the infinitude of prime numbers 
which was shown by Euclid, Sieve of Eratosthenes for prime numbers, consideration of 
Diophantine equations) led to the birth of elementary theory of numbers. Later investiga- 
tions (e.g. Fermat's Last theorem, investigations of zeta-function by Euler and Riemann) 
have led to the development the great branch of mathematics - the analytical theory of 
numbers. The important problems of modern number theory are the problems of effec- 
tive search of prime numbers and complexity (and factorization) or primality of arbitrary 
natural number. Today essential efforts are applied for consideration of these questions: 
many fast methods for factorization (e.g. Pollard p-algorithm [I], elliptic curve method 
[2], quadratic sieve algorithm [3], general number field sieve [4]) and tests of primality 
(e.g. Miller- Rabin test [5], Lucas-Lehmer test [6], AKS test [7]) are developed, projects 
on search of large prime numbers are running (GIMP (8], PrimeGride [9], Seventeen or 
Bust [10] and others), databases of prime numbers are updated (e.g. The List of Largest 
Known Primes [11], OEIS [12]). Such interest is due to not only desire to understand 
the nature of numbers and to resolve some fundamental problems of mathematics but 
also due to opportunity of practical application of properties of natural numbers, e.g. 
for the cryptographic systems with public key (for example RSA [13]), hash-tables, and 
generation of pseudorandom numbers (e.g. Mersenne twister [14]). 

The aim of this paper is to suggest a theoretical instrument for investigation of 
quantity of divisors of natural number, which may be useful for further investigations of 
properties of natural numbers. 



2 Integral representation of divisor function 

Main idea is similar to the one of method of trigonometric sums [15] . Using the fact that 
exp (2nix) equals unity only for integer value of real number x one may say that function: 

f s (x) = exp (2nix) + exp {2nis/x) — 2 (1) 

has zeros only for values of x which are divisors of natural number s (we will denote these 
numbers as Xj, where superscript s highlights dependence on natural number s, subscript 
j denotes integer sequence number of corresponding zero). Note that 

exp (2irix s j ) = 1, exp (27ris/Xj) = 1 (2) 

Unfortunately, direct use oi f s (x) (e.g. search of global minima of f s {x) by using numerical 
iteration methods) for factorization of natural number s or for counting number of divisors 
of s is almost impossible because of frequent oscillations of function f s (x) (see fig. [lj. 

Number of zeros of this function in the interval [1; s] gives us the number of divisors 
of s. To count the number of zeros of f s one can use logarithmic residue of f s . In general 
case the logarithmic residue is the following (e.g. [16J): 



where / is a function of complex number z, L is a closed contour on the complex plane 
C, Nf is the number of zeros of function / within the contour L (taking into account the 
orders of zeros), Pf is the number of poles of function / within the contour L (taking into 
account the orders of poles). If function / has no poles within contour L (e.g. / is the 
smooth functionof z) then one may write: 

7 [/,£]= Yl k * ( 4 ) 

Zj 6L:/( Zj )=0 

where L is the field within the contour L (i.e. L is the field limited by the contour L on 
the complex plane C), Zj are zeros of function / in the field L, kj is the order of zero Zj, 
i.e. in the vicinity of Zj the function / (z) can be presented as 

/ (z) = ^-T^ (* ~ z rf j + ((* ~ *i) kj+1 ) (5) 

Using flU) one can construct the expression for counting number of zeros of function 
f a (z) in the interval [1; y/s\. Let us define the following functions: 

r_( a ) = 7 [/.,^r] , r c (s) = 7 [f s , L°] , r + (s) = 7 [/., L+] (6) 
where L~ is the contour given (e.g.) by the following points z (see fig. EJ): 

z = l-e + iy, ye [-e; +e] 
x + ie, x G [(1 - e) ; (y/s - e)] 

(a/s -e)+iy, ye [+e; -e] 
x — ie, a; 6 [(\/s — e) ; (1 — e)] (7) 

is the contour given (e.g.) by the following points z (see fig. |2|: 



z = y/s - e + zy, ye [-e; +ej 
x + ie, x e [(y/s — e) ; (-y/s + e)] 
(v 7 ^ + e) + zj/, y G [+e; -e] 
x-ie, x 6 [(Vs + e) ; (Vs - e)] (8) 

is the contour given by the following points z (see fig. |2J): 

z = l-e + iy, ye [-e; +e] 
x + ze, x G [(1 - e) ; + e)] 

(a/s + e) + zj/, j/ G [+e; -e] 
x-z'e, xG [(Vs + e);(l-e)] (9) 

where e is the positive real number. Parameter e should be small enough to exclude zeros 
(z^ denote zeros of f s (z) in complex plane) which do not lie on the real axis. Note that 
r+(s) = r_(s) + r c (s). 
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Studying values of t_(s), t c (s) and r + (s) one should calculate the orders of zeros 

.r 



Zj which lie on the real axis (i.e. Zj 



fs(z): 



f's (*) = 2m 



exp (2mz) exp 



This demands calculation of derivative of 
2ms' 



(10) 



Taking into account (j2J) one may write 



f's 



2m 



Thus f' s \Xj) f for x s j G [1; \/s), and therefore all zeros Xj in interval [1; yfs) have order 
which equals unity. If yfs is also divisor of s then one should investigate the order of this 
zero. From ( JTT1) it is evident that this order is not equal unity. Thus one should find the 
second order of f s (z): 



fs CO = 2™ 



2s ( 2iris 
2ki exp \2kiz) H — - exp 



s ( 2ms 
+ 2ki— exp 

\ 2; 



12] 



If the value yfs is an integer number (this case is interesting for us) then one may write: 

f'J (^i) = -8tt 2 + 47TZS" 1 / 2 ^ (13) 



Expression ( 1T3|) means that the order of zero in point y/s equals 2 (if yfs is integer and 
divisor of number s) or (if y/s is not integer and not a zero of function f s (z)), i.e. 
r c (s) = or 2. 

Knowledge of behaviour of functions T_(s), t c (s) and r + (s) allows us to calculate 
the number of divisors r (s) of integer number s. Taking into account symmetry of 
quantity of divisors of s relative to the number yfs one may write: 



t(s) = 2t-(s) + -t c (s) 



2rJs) 



-TAS 



Explicit expressions for r_(s), t c (s) and t+(s) are the following: 

exp {2Kiz) — (s/z 2 ) exp [2Kisj z) 



T-(S 



T+(S) 



L - exp (2mz) + exp (2ms / z) — 2 

exp(2mz) — (s/z 2 ) exp (2ms /z) 
L c exp (2mz) + exp (2Kis/z) — 2 

exp (2Kiz) — (s/z 2 ) exp (2kis/ z 



L + exp (2mz) + exp (2ms/ z) — 2 



-dz 



dz 



dz 



(14) 

(15) 

(16) 
(17) 



3 Prime number equation 

One of the interesting cases of using representation (I14II17P appears when s is a prime 
number. In this case t(s) = 2 and one can write the equation for prime numbers using 
equivalent condition t + (s) = 1. To prove the equivalence of r(s) = 2 and r + (s) = 1 let 
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us consider the following: 

1) From t(s) > 2 and r(s) < 2r + (s) (from (JHJ)) one can show that r+(s) = 1 =>• r(s) = 2, 

2) From r+(s) > 1, T+(s) 7^ 2, and r + (s) < r(s) one can show that r(s) = 2 =>- 7+(s) = 1. 

Thus one can write the following equation for prime numbers p: 

f exp (2mz) - (p/z 2 ) exp (2mp/z) ^ = ^ 
J L + exp (27riz) + exp (2irip/z) — 2 



4 Conclusion 

Formulae (I14H17P for an integral representation of divisor function and equation ([15]) for 
prime numbers have been shown in the present paper. Unfortunately the calculation of 
expressions (I14H17P and ([15]) for the typical values of s and p (which are needed in practical 
problems) demands huge computational time and resources. Thus direct application of 
formulae ( 11411171) and (I18jl is impossible. Nevertheless, these formulae may be useful 
for further investigations of properties of natural numbers and solving of practical and 
fundamental problems of number theory such as construction of new tests of primality 
and proof of the infinitude of twin primes. 
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Figure 1: Dependence of module of function f s (x) on real variable x. 
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Figure 2: Examples of contours L in complex plane for calculation of logarithmic residue 
of function f s (x). 
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